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BACKSCAT  Lidar  Backscatter  Simulation 
User’s  Manual  for  Version  1.0 


1  INTRODUCTION 

Backscatter  from  atmospheric  aerosols  can  produce  significant  returns  in  monostatic  laser 
radar  (lidar)  systems.  The  Air  Force  Geophysics  Laboratory  (AFGL)  is  developing  lidar  systems 
that  will  measure  this  backscatter.  To  aid  in  the  design  and  use  of  such  systems,  SPARTA  has 
developed  a  simulation  program,  BACKSCAT,  that  calculates  the  backscatter  return  for  various 
lidar  systems,  viewing  aspects,  and  atmospheric  conditions.  The  preliminary  version  that  is 
described  in  this  report  was  designed  to  look  at  aerosol  effects. 

The  BACKSCAT  computer  simulation  code  is  designed  to  be  easy  to  use  while  providing  the 
operator  with  maximum  flexibility.  All  simulation  and  configuration  parameters  can  be  changed 
interactively  through  an  integrated  menu  interface.  The  simulation  output  is  written  in  tabular 
form  to  a  data  file  that  can  be  read  by  standard  plotting  programs.  A  log  file  is  also  generated 
during  the  execution  of  the  simulation  that  provides  for  verification  of  simulation  inputs  and 
traceability  of  results. 
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I.l  Organization  of  Report 

This  manual  is  a  user’s  guide  to  the  operation  of  the  BACKSCAT  lidar  backscatter  simulation 
code.  Section  2  summarizes  key  features  of  the  code  itself.  Section  3  describes  the  steps  required 
to  install  the  code  on  a  computer  system.  Section  4  discusses  the  interactive  execution  of  the 
code.  Section  5  provides  samples  of  the  output  generated  by  the  simulation.  Finally,  Section  6 
provides  recommendations  for  future  modifications. 

Three  appendices  contain  supplementary  information  about  BACKSCAT.  Appendix  A  sum¬ 
marizes  the  program  modules  that  make  up  the  code.  Appendix  B  discusses  the  adaptations 
which  may  be  necessary  when  installing  the  code  on  computers  with  different  operating  sys¬ 
tems.  Finally,  Appendix  C  describes  the  format  and  syntax  of  the  files  used  by  BACKSCAT. 
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2  KEY  FEATURES  OF  BACKSCAT 
2.1  Computer  System  Portability 

This  version  of  BACKSCAT  was  developed  on  an  IBM-compatible  personal  computer.  The 
code  was  written  and  compiled  using  Microsoft™  Fortran  version  4.01  of  Fortran  77.  Being 
written  in  Fortran  77,  BACKSCAT  should  be  able  to  run  on  virtually  any  host  computer.  Minor 
changes  may  be  necessary  when  installing  the  code  on  other  computer  systems  due  to  differences 
in  operating  systems,  but  every  reasonable  effort  has  been  made  to  isolate  these  changes  to  a 
minimal  number  of  subroutines.  Appendix  A  contains  a  list  of  all  of  the  programs  that  make  up 
BACKSCAT  and  a  brief  description  of  what  they  do. 


2.2  Summary  of  the  Physics 

BACKSCAT  simulates  the  return  from  a  monostatic  lidar  through  a  solution  of  the  lidar 
equation  (Cf.  Klett1).  The  backscatter  return  observed  by  a  lidar  system  operating  at  the 
wavelength  A  is  given  by 


P*(A,r)  =  +  /MA,r 


-2/[cra(A,r)  +  am(A,r)]dr 
0 


where: 

PrM  is  the  power  received  from  range  r, 
t  is  the  pulse  length, 
c  is  the  speed  of  light, 

Ar  is  the  collecting  area  of  the  receiver  aperture, 

Pj  is  the  transmitted  power, 

0a{ A,r)  is  the  aerosol  backscatter  coefficient  at  range  r, 

Pm{ A, r)  is  the  molecular  backscatter  coefficient  at  range  r, 
aa(A,r)  is  the  aerosol  attenuation  coefficient  at  range  r,  and 
am(A,r)  is  the  molecular  attenuation  coefficient  at  range  r. 

The  lidar  return  can  be  normalized  by  dividing  by  the  transmitter  power  to  obtain  the 
fractional  power  returned  as  a  function  of  range.  The  lidar  return  can  also  be  range  compensated 
by  multiplying  by  the  square  of  the  range  to  obtain  a  measure  of  the  actual  aerosol  levels. 

The  simulation  can  be  run  with  molecular  backscatter  and  attenuation  contributions  included 
or  excluded.  One  can  also  selectively  include  molecular  backscattering  and  neglect  molecular 
absorption  and  vice  versa. 


1  Klett,  J.  D.  (1981)  Stable  Analytical  Inversion  Solution  for  Processing  Lidar  Returns,  Applied 
Optics  20:211. 
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BACKSCAT  contain*  a  database  of  attenuation  and  backscatter  coefficients  used  in  the 
AFOL  aerosol  models2.  The  user  can  also  supply  a  specific  set  of  aerosol  and  molecular 
attenuation  and  backscattering  coefficients.  The  BACKSCAT  simulation  obtains  the  attenuation 
and  backscatter  coefficients  by  linear  interpolation  of  an  altitude  index  into  a  propagation  profile 
table.  The  altitude  at  a  particular  range  is  obtained  from  the  “flat  earth’ ’  relationship 

h{r)  -  h0  +  rsln(^) 

where: 

h(r)  is  the  altitude  at  range  r, 
ho  is  the  sensor  altitude,  and 
^  is  the  sensor  viewing  elevation  angle. 

The  area  of  the  receiver  aperture.  Ap,  is  computed  from  the  specified  aperture  and  obscu¬ 
ration  diameters,  assuming  a  standard  circular  aperture  with  a  single  circular  obscuration. 


2  R.W.  Fenn  et  al.  (1985),  “Optical  Infrared  Properties  of  the  Atmosphere”,  Chapter  18  in 
Handbook  of  Geophysics  and  the  Space  Environment ,  A.S.  Jursa  (ed.),  Air  Force  Geophysics 
Laboratory,  Hans  com  Air  Force  Base,  Massachusetts,  AFGL-TR-85-0315,  ADA167000,  and 
E.P.  Shettle  (1988)  Private  Communication. 
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3  INSTALLATION 

3.1  Installation  on  IBM  or  Compatible  Computers 

BACKSCAT  is  designed  to  operate  on  an  IBM  or  compatible  personal  computer  with  an 
8087  math  co-processor  and  operating  with  PC-DOS  or  MS-DOS.  The  code  can  be  run  on 
computers  with  dual  floppy  disk  drives  or  a  hard  disk.  The  use  of  hard  disk  is  recommended. 

To  install  the  code,  insert  the  diskette  containing  the  executable  program  into  a  5$  inch 
diskette  drive  and  then  type  the  command: 

[dl :  ]  INSTALL  dl :  [<42:1  [directory) 

and  then  press  the  ENTER  or  RETURN  key.  (For  simplicity  sake,  we  shall  only  make  reference 
to  the  ENTER  key  in  this  report  but  its  use  is  the  same  as  the  RETURN  key.)  In  the  above 
command,  dl :  is  the  diskette  drive  containing  the  programs,  d2 :  is  the  destination  drive,  and 
directory  is  the  directory  or  subdirectory  into  which  the  program  is  to  be  installed.  (The 
components  of  the  INSTALL  command  line  shown  in  brackets  are  optional  and  the  brackets  are 
not  typed  when  the  command  is  entered).  If  a  destination  directory  is  not  specified,  the  current 
default  directory  will  be  assumed.  If  a  specific  destination  directory  is  desired,  it  must  be  created 
prior  to  installation.  (If  the  program  is  to  be  customized,  the  source  (.FOR)  files  must  also  be 
copied  from  the  source  distribution  diskette  using  die  PC-DOS  or  MS-DOS  COPY  command.) 

If  an  error  is  made  during  the  installation  procedure,  such  as  from  neglecting  to  add  the  dl : 
after  the  command  INSTALL,  an  error  message  will  be  displayed  on  the  monitor.  Correct  the 
installation  commands  that  were  entered  and  press  the  ENTER  command. 

3J  Installation  on  Non  PC  or  MS-DOS  Computers 

Installation  of  the  BACKSCAT  program  on  a  computer  using  another  operating  system  is 
slightly  more  complex.  The  source  and  program  data  files  (all  files  on  the  source  distribution 
diskette)  must  first  be  transferred  from  the  source  program  distribution  diskette  to  the  host 
computer  system.  The  source  files  must  then  be  compiled  and  linked  on  the  host  system  to  create 
the  correct  executable  code.  It  may  additionally  be  necessary  to  rename  files  and/or  modify  some 
subroutines  to  accommodate  operating  system  differences,  as  described  in  Appendix  B. 

The  DOS  file  ANSII.SYS  file  should  be  added  to  the  AUTOEXEC.BAT  file  for  the  computer 
on  which  BACKSCAT  will  be  run.  This  will  enable  the  code  to  take  full  advantage  of  the  menu 
operations  of  the  ptogram. 
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4  PROGRAM  OPERATION 
4.1  Beginning  the  Program 

The  BACKSCAT  program  may  be  executed  under  PC-DOS  or  MS-DOS  by  entering  the 
program  name  (BACKSCAT)  at  the  DOS  system  prompt  and  then  press  ENTER.  The  required 
commands  for  executing  BACKSCAT  on  computers  with  different  operating  systems  may  differ 
and  the  user  is  urged  to  consult  the  appropriate  operating  system  manual  for  specific  instructions. 

BACKSCAT  first  displays  a  title  screen,  as  shown  in  Figure  1,  with  a  prompt  for  a  system 
configuration.  The  user  must  enter  the  name  of  an  existing  configuration.  A  configuration 
named  DEFAULT  is  supplied  with  the  program  and  must  be  entered  the  first  time  BACKSCAT 
is  used.  Additional  configurations  may  be  created  by  the  user  within  die  program  as  described 
below  and  subsequently  used  for  program  initialization.  The  program  converts  the  configuration 
name  to  a  file  name  and  attempts  to  read  the  information  in  the  file.  If  the  program  cannot 
locate  the  specified  configuration  file,  a  message  will  be  displayed  prompting  the  user  to  enter 
another  configuration  name.  The  user  is  cautioned  not  to  change  the  conditions  specified  in  the 
DEFAULT  configuration  file. 

The  configuration  file  name  also  contains  names  of  filer,  that  describe  the  lidar  system  and 
atmospheric  structure  to  be  modeled.  If  either  file  cannot  be  located,  the  user  is  prompted  for 
another  name. 

4 1  Menu  Control  System 

BACKSCAT  utilizes  a  series  of  master  menus  to  define  a  set  of  simulation  parameters. 
These  parameters  define  die  lidar  system,  atmospheric  structure,  and  configuration  conditions. 
These  master  menus,  in  turn,  link  to  a  series  of  submenus  that  permit  the  user  to  define  the 
individual  parameters. 

Each  master  menu  or  submenu  is  divided  into  a  status  area  and  a  dialog  area  separated  by  a 
row  of  asterisks.  The  status  area  displays  the  current  set  of  conditions.  The  dialog  area  contains 
prompts  to  aid  the  user  in  changing  any  of  the  operating  conditions. 

4 3  Main  Program  Status 

Once  the  desired  configuration  has  been  loaded,  the  program  will  display  a  Main  Status 
Menu,  as  shown  in  Figure  2,  showing  the  conditions  of  the  specified  configuration.  The  values 
shown  in  the  figure  and  subsequent  figures  correspond  to  the  default  configuration. 

The  Main  Status  Menu  displays  the  names  of  the  files  containing  the  information  about  the 
lidar  system,  the  atmospheric  model,  the  propagation  profile,  the  molecular  absorption  profile, 
the  simulation  log,  the  simulation  output,  and  the  mode  of  propagation.  The  dialog  area  at  the 
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Figure  1.  Opening  Menu  of  BACKSCAT  Displaying  the  Prompt  for  the  Name  of  the  Desired 
Configuration 

bottom  of  the  menu  contains  letters  that  one  can  use  to  change  any  of  the  above  parameters.  If 
one  does  not  enter  a  given  letter,  the  program  will  automatically  scroll  through  each  parameter. 
The  monitor  will  display  the  current  values  for  the  given  parameter  and  a  prompt  to  enter  the 
new  information.  Text  entries  can  be  entered  in  upper  or  lower  case  letters.  The  program  will 
automatically  convert  the  text  to  upper  case.  If  one  does  not  wish  to  change  a  given  parameter, 
press  the  ENTER  key. 


4A  Quitting  Menus  and  Program  Termination 

The  “Q  -  Quit”  option  on  any  menu  universally  passes  control  out  of  that  menu.  From 
any  menu  other  than  the  Main  Status  Menu,  control  reverts  to  the  next  higher  menu.  When  the 
“Q  --  Quit”  option  is  selected  from  the  Main  Status  Menu,  the  program  terminates  and  returns 
control  to  the  operating  system.  Menus  that  do  not  include  a  "Q  —  Quit”  option  (generally 
used  where  the  user  must  select  a  value  from  a  discrete  set  of  valid  responses)  will  automatically 
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Figure  2.  Main  Status  Menu  Displaying  the  Conditions  of  the  Specified  Configuration 

return  control  to  the  next  higher  menu  when  the  selection  is  completed. 

When  the  user  chooses  to  quit  any  of  the  master  menus,  BACKSCAT  first  checks  to  see  if  any 
changes  were  made  to  the  affected  data  set  (lidar  system,  atmospheric  structure,  or  configuration). 
If  changes  were  made,  the  quit  submenu  is  displayed  an4  the  user  may  choose  ”C  -  Cancel" 
to  cancel  the  changes,  ”K  --  Keep”  to  keep  the  changes  for  the  current  session  without  saving 
to  a  file,  or  ”S  -  Save”  to  save  the  modified  data  set  to  a  file.  If  the  user  chooses  to  save 
the  modified  data  set,  the  current  name  for  the  data  sets  (lidar  system,  atmospheric  structure,  or 
configuration)  is  used,  overwriting  any  existing  data  set  with  the  same  name. 


DEFAULT 
TUT 
(MOHS) 
DEFAULT 
DEFAULT 
ATMOS.  MODEL  /  ON 
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Figure  3.  The  Master  Lidar  System  Menu 
4J  Changing  Lidar  Systems 

The  parameters  governing  the  lidar  system  may  be  changed  by  selecting  the  “L  --  Lidar 
System”  option  from  the  Main  Status  Menu.  This  option  will  result  in  the  display  of  the  Master 
Lidar  System  Menu  as  shown  in  Figure  3.  To  change  a  given  parameter,  either  type  the  letter 
corresponding  to  the  parameter  and  press  ENTER  or  press  ENTER  until  the  program  scrolls  to 
the  desired  parameter. 

4.5.1  ”N  -  Name  of  System”  Option 

This  option  prompts  the  user  for  a  new  lidar  system  name.  The  current  name  will  be 
displayed  along  with  a  prompt  for  the  new  name.  Type  in  the  new  name,  using  eight  characters 
or  less,  and  press  ENTER.  If  you  do  not  wish  to  make  a  change,  merely  press  ENTER.  The 
name  entered  will  be  displayed  on  the  Master  Lidar  System  Menu. 


N  --  Naae  of  Syatea  L  —  Load  Syatm  T  —  Traneaitter 

*  Receiver  •  —  Signal  Proceeeor  H  --  sanaor  Haight 
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4.5.2  ”L  ~  Load"  Option 

This  option  converts  the  lidar  system  name  to  a  file  name.  If  the  file  exists,  the  program  then 
loads  the  lidar  system  information  contained  in  the  file.  If  the  file  does  not  exist,  the  information 
relating  to  the  current  lidar  system  is  retained  and  referenced  by  the  new  name. 

4.5.3  *‘T  -  Transmitter”  Option 

This  option  displays  the  Lidar  Transmitter  Submenu,  as  shown  in  Figure  4,  displaying  the 
current  transmitter  parameters.  The  parameters  that  can  be  changed  are  the  lidar  wavelength, 
in  microns;  the  pulse  energy,  in  joules;  and  the  pulse  duration,  in  microseconds.  To  change  a 
transmitter  parameter,  type  the  letter  of  the  parameter  and  press  ENTER  or  press  ENTER  until 
the  desired  parameter  is  displayed  in  the  dialog  area.  To  return  to  the  Master  Lidar  System 
Menu,  type  ”Q  --  Quit”  and  press  ENTER. 


Figure  4.  The  Lidar  Transmitter  Submenu 
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4.5.4  "R  -  Receiver"  Option 

This  option  displays  the  Lidar  Receiver  Submenu,  as  shown  in  Figure  5,  showing  the  current 
receiver  parameters.  The  parameters  that  can  be  changed  are  the  receiver  aperture  diameter, 
in  centimeters,  and  the  receiver  obscuration  diameter,  in  centimeters.  To  change  one  of  the 
parameters,  type  the  letter  of  the  choice  and  press  ENTER  or  press  ENTER  until  the  parameter 
appears  in  the  dialog  area.  To  return  to  the  Master  Lidar  System  Menu,  type  Q  and  press 
ENTER. 


r 

Aperture  Disaster  (Centlaetere) 

100.0 

Obscuration  Disaster  (Centlaetere) 

0.0 

RECEIVER : 

0  --  Obscuration  Disaster 

A  *•  Aperture  Disaster 

0  —  Quit 

^  Enter  Selection,  or  Press  ENTER  Ray  for  Default 

(Hifhllqt.ted)  Selection. 

Figure  5.  The  Lidar  Receiver  Submenu 


4.5.5  "S  -  Signal  Processor"  Option 

This  option  displays  the  Lidar  Signal  Processor  Submenu,  as  shown  in  Figure  6,  and  the 
current  signal  processor  parameters  menu.  The  signal  processor  parameters  are  the  nearest  or 
minimum  range  that  can  be  detected,  in  kilometers,  and  the  farthest  or  maximum  range,  in 
kilometers.  To  change  a  parameter,  enter  the  letter  displayed  and  press  ENTER  or  press  ENTER 
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Figure  6.  The  Udar  Signal  Processor  Submenu 

until  the  parameter  appears  in  the  dialog  area.  To  return  to  the  Main  Lidar  System  Menu,  type 
Q  and  press  ENTER. 

4.5.6  Sensor  Location  and  Viewing  Aspect 

The  user  may  select  “H  -  Senior  Height”  to  specify  the  altitude  of  the  sensor  above  the 
ground.  Similarly,  the  user  may  specify  the  viewing  direction  by  selecting  ‘‘A  -  Azimuth 
Angle”  or  ‘‘E  -•  Elevation  Angle”  to  specify  the  azimuth  and  elevation  angles,  respectively. 

The  altitude  must  be  entered  in  kilometers,  and  both  angles  must  be  entered  In  degrees. 
The  elevation  angle  is  measured  in  degrees  above  (positive)  or  below  (negative)  the  horizontal 
plane.  The  azimuth  angle  currently  has  no  effect,  but  is  included  for  use  in  future  upgrades  to 
the  simulation. 
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4.5.7  ”Q  -  Quit”  Option 

This  option  exits  the  Master  Lidar  System  Menu  and  returns  control  to  the  Main  Status 
Menu.  Type  Q  and  press  ENTER. 

4.6  Atmospheric  Model 

The  atmospheric  structure  specifications  for  the  AFGL  atmospheric  model  may  be  changed 
by  selecting  the  “A  -  Atmospheric  Model”  option  from  the  Main  Status  Menu.  This  option 
displays  the  Master  Atmospheric  Model  Menu,  as  shown  in  Figure  7. 

( - " 


Name  of  Structure 

DEFAULT 

Upper  Atsosphere  --  Aeroeol  Level 
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Figure  7.  Master  Atmospheric  Model  Menu 


In  the  AFOL  atmospheric  models,  the  atmosphere  is  divided  into  four  layers.  From  the 
ground  upward,  these  layers  are  the  boundary  layer,  the  troposphere,  the  stratosphere,  and  the 
upper  atmosphere.  The  heights  of  each  of  the  layer  boundaries  are  currently  preset  but  can  be 
set  by  the  user  if  an  appropriate  data  base  is  supplied. 
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The  data  base  that  is  required  is  a  table  of  extinction  coefficients  for  a  wavelength  of  O.SS  n m 
as  a  function  of  altitude  for  one  of  the  AFGL  reference  aerosol  models.  This  data  base  contains 
the  information  on  how  to  scale  the  aerosol  abundances  as  a  function  of  altitude.  The  default 
values  that  are  provided  with  BACKSCAT  are  contained  in  a  file  called  STANDARD.SCL.  The 
data  in  the  file  are  read  in  with  the  FORTRAN  format  (F6.1.9E8.2).  Each  record  of  the  file 
contains  the  altitude,  in  km,  and  the  extinction  coefficients  to  different  conditions.  (The  precise 
columnar  format  of  the  data  file  are  based  on  Table  18- 10a.  from  Fenn  et  alA) 

4.6.1  “N  -  Name  of  Structure”  Option 

This  option  prompts  the  user  for  a  new  atmospheric  structure  name.  To  enter  a  new  name, 
type  a  name  of  eight  characters  or  teas,  and  press  ENTER.  The  name  entered  becomes  the  name 
of  the  current  atmospheric  structure. 

4.6.2  ‘‘L  -  Load”  Option 

This  option  converts  the  atmospheric  structure  name  to  a  file  name.  If  the  file  exists,  the 
program  then  loads  the  information  about  the  atmospheric  structure  contained  in  the  file.  If  the 
file  does  not  exist,  the  information  about  the  currently  displayed  atmospheric  structure  is  retained 
and  referenced  by  the  new  name. 

4.6.3  *‘B  -  Boundary  Layer”  Option 

This  option  is  used  to  describe  the  aerosol  conditions  in  the  planetary  boundary  layer.  The 
conditions  under  user  control  are  the  type  of  aerosol,  the  surface  value  of  the  relative  humidity, 
and  the  surface  value  of  visibility.  This  option  displays  the  boundary  layer  submenu.  To  select 
this  option,  type  B  and  press  ENTER. 

4  6.3.1  “A  ACKBflr  Option 

This  option  allows  the  user  to  describe  the  type  of  aerosols  in  the  boundary  layer.  The  user 
can  select  ‘‘R  -  Rural”  for  rural  aerosols,  ”U  -  Urban”  to  urban  aerosols,  ”M  --  Maritime”  to 
maritime  aerosols,  “O  -  Oceanic”  to  oceanic  aerosols,  ”T  -  Tropospheric”  to  tropospheric 
aerosols,  or  ”F  ~  Fog”  to  fog  conditions.  Type  in  the  appropriate  letter  and  press  ENTER. 

If  the  user  selects  fog  conditions,  the  fog  submenu  is  displayed.  The  user  can  then  select 
"A  -  Advecdve”  for  advective  fog  or  ”R  --  Radiative”  to  radiative  fog.  Type  A  or  R  and  then 
press  ENTER.  Upon  completion  of  the  selection,  control  returns  to  the  boundary  layer  aerosol 
submenu. 
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4.6.3.2  ”fl~  Humidittl’.Qprion 

This  selection  allows  the  user  to  specify  the  relative  humidity  at  the  surface  of  the  boundary 
layer.  Type  in  the  desired  value,  in  percent  (%),  and  press  ENTER. 

4.6.3.3  “V  -«  Visibility”  Potion 

This  selection  prompts  the  user  for  the  visibility  at  the  surface  of  the  boundary  layer,  which 
is  then  used  to  calculate  the  extinction  coefficient  at  the  surface.  To  enter  the  visibility,  enter 
a  positive  value,  in  km.  One  can  also  assume  that  the  extinction  coefficient  in  the  boundary 
layer  is  uniform  throughout  the  boundary  layer.  To  invoke  this  option,  enter  a  negative  value, 
in  km-1.  Press  ENTER  after  either  value  has  been  typed  in. 

4.6.3.4  “0  -  Quit**  Potion 

This  selection  returns  control  to  the  Master  Atmospheric  Model  Menu.  Type  Q  and  press 
ENTER. 

4.6.4  “T  Tropospheric  Humidity”  Option 

This  option  prompts  the  user  for  the  relative  humidity  of  the  troposphere.  Type  in  the  desired 
value,  in  percent  (%),  and  press  ENTER. 

4.6.5  ‘‘S  -*  Stratosphere”  Option 

The  Stratosphere  Option  is  used  to  describe  the  aerosol  conditions  in  the  stratosphere.  The 
conditions  that  must  be  specified  are  the  type  and  profile  of  the  stratospheric  aerosols.  To  select 
this  option,  type  S  and  press  ENTER. 

46.5.1  **E  Extinction"  Option 

This  selection  displays  the  choices  for  the  stratospheric  extinction.  The  user  can  choose  ”B  - 

-  Background”  for  background  extinction  conditions,  “F  -  Fresh  Volcanic”  for  fresh  volcanic 
extinction  conditions,  or  “A  -  Aged  Volcanic”  for  aged  volcanic  extinction  conditions.  Type 
B,  F  or  A  and  press  ENTER.  Control  will  then  return  to  the  stratospheric  submenu. 

4.6.5.2  ”P  ~  Profile”  Potion 

This  selection  displays  the  choices  for  the  stratospheric  profile.  The  user  can  choose  ”B  - 

-  Background”  for  a  background  stratospheric  aerosol  profile,  ”M  -  Moderate  Volcanic”  for 
a  moderate  volcanic  aerosol  profile,  ”H  ••  High  Volcanic”  for  a  high  volcanic  aerosol  profile, 
or  ”E  -  Extreme  Volcanic”  for  an  extreme  volcanic  aerosol  profile.  Type  in  the  letter  of  the 
desired  choice  and  press  ENTER.  Control  returns  to  the  stratospheric  menu  after  a  selection  is 
made. 
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4.6.5. 3  “0  --Qutr-Qpaon 

This  selection  returns  control  to  the  Master  Atmospheric  Model  Menu.  Type  Q  and  press 
ENTER. 

4.6.6  “U  --  Upper  Atmosphere”  Option 

This  selection  displays  the  options  for  the  upper  atmosphere  aerosols.  The  user  can  choose 
”N  -  Normal”  for  normal  upper  atmospheric  aerosol  conditions  or  ”E  ~  Extreme”  for  extreme 
meteoric  aerosol  conditions.  Type  N  or  E  and  then  press  ENTER.  After  a  selection  is  made, 
control  returns  to  the  Master  Atmospheric  Model  Menu. 

4.6.7  ”D  -•  Distribution”  Option 

This  selection  displays  the  choices  for  the  seasonal  aerosol  distributions.  The  user  can  select 
from  a  Fall/Winter  aerosol  profile  or  Spring/Summer  aerosol  profile.  Type  F  or  S  and  then  press 
ENTER.  Control  then  returns  to  the  Master  Atmospheric  Model  Menu. 

4.f.8  ”H  --  Heights  of  Layers”  Option 

This  option  displays  the  Aerosol  Layer  Height  Submenu,  as  shown  in  Figure  8.  and  allows 
the  user  to  change  the  boundary  heights  between  different  types  of  aerosols.  The  user  may  select 
*‘B  -Boundary  Layer”  to  set  the  upper  limit  of  the  boundary  layer,  “T  -  Troposphere”  to  set 
the  upper  limit  of  the  troposphere,  ”S  --  Stratosphere”  to  set  the  upper  limit  of  the  stratosphere 
”U  -*  Upper  Atmosphere”  to  set  the  upper  limit  of  the  assumed  upper  atmosphere,  or  “Q  - 
Quit”  to  return  control  to  the  Master  Atmospheric  Model  Menu.  Type  in  the  letter  corresponding 
to  one’s  choice  and  press  ENTER.  The  program  will  prompt  the  user  to  enter  a  new  boundary 
layer  height  in  kilometers.  The  lower  limit  of  the  boundary  layer,  of  course,  is  zero.  The  lower 
limit  of  the  other  layers  is  the  upper  limit  of  the  next  lower  layer. 

The  values  for  layer  heights  that  are  supplied  in  the  default  configuration  correspond  to  those 
in  the  AFGL  model  atmospheres  and  correspond  to  the  scaling  factors  discussed  in  Section  4.6. 
If  these  heights  are  to  be  changed,  the  user  must  supply  a  new  set  of  aerosol  scale  factors  which 
are  contain  the  new  ret  of  layer  heights.  Or,  the  code  must  be  modified  to  include  a  procedure 
for  interpolating  the  scale  factors  for  the  user  supplied  layer  heights. 

4.6.9  “Q -- Quit”  Option 

This  option  exits  the  Master  Atmospheric  Model  Menu  and  returns  control  to  the  Main  Status 
Menu.  Type  Q  and  then  press  ENTER. 
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Figure  8.  The  Aerosol  Layer  Height  Submenu 
4.7  Propagation  Profile  Name 

If  “P  --  Propagation  Profile"  is  selected  from  the  Main  Status  Menu,  the  user  is  prompted 
for  a  new  propagation  profile  name.  Type  in  a  name  of  eight  characters  or  less  and  press  ENTER. 

When  the  simulation  is  executed  from  the  Main  Status  Menu,  as  described  below,  the  prop¬ 
agation  profile  name  is  convened  to  a  file  name.  If  the  propagation  mode  is  set  for  propagation 
profile,  the  corresponding  file  is  read  to  obtain  the  propagation  data  and  the  atmospheric  model  is 
ignored.  If  the  propagation  mode  is  set  for  atmospheric  model,  a  propagation  profile  is  generated 
by  the  atmospheric  model  for  the  current  atmospheric  structure  and  written  to  the  corresponding 
file,  overwriting  any  existing  file  of  the  same  name. 

When  the  propagation  mode  is  selected  to  the  propagation,  so  that  the  simulation  will  read 
a  propagation  profile  from  a  file,  the  file  must  be  written  in  standard  tabular  form.  The  columns 
of  the  propagation  profile  must  be:  (1)  altitude  (kilometers),  (2)  aerosol  extinction  coefficient 
(inverse  kilometers),  (3)  aerosol  absorption  coefficient  (inverse  kilometers),  (4)  aerosol  scatter 
coefficient  (inverse  kilometers),  (5)  aerosol  backscatter  coefficient  (inverse  meter-steradians),  (6) 
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molecular  extinction  coefficient  (inverse  kilometers),  and  (7)  molecular  backscatter  coefficient 
(inverse  meter-steradians). 

The  term  “standard  tabular  form”  means  that  each  row  of  data  correspond  to  a  separate 
altitude.  The  data  are  read  in  “free  format”  form.  Altitudes  must  appear  in  increasing  order. 
Although  there  is  no  limit  to  the  number  of  altitudes  that  may  appear  in  a  data  file,  the  current 
version  of  BACKSCAT  will  only  read  the  first  one  hundred  lines  of  data. 

4.8  Molecular  Absorption 

The  “M  •  Molecular  Absorption”  option  allows  the  user  to  specify  a  molecular  absorption 
profile  for  use  by  the  lidar  equation.  If  the  user  does  not  wish  to  include  molecular  absorption, 
type  in  (NONE)  or  NONE.  If  the  user  wishes  to  include  molecular  absorption,  enter  in  the  name 
(without  any  file  extensions)  of  the  file  containing  the  absorption  data. 

The  Ale  containing  the  molecular  absorption  must  contain  two  columns  of  data;  the  first 
column  is  the  altitude  in  kilometers  and  the  second  column  is  the  absorption  coefficient  in 
inverse  kilometers.  The  data  must  be  arranged  such  that  the  alutudes  are  in  increasing  order. 

4.9  Simulation 

The  simulation  option  may  be  entered  by  selecting  “S  -  Simulation”  from  the  Main  Status 
Menu.  The  user  will  then  be  given  the  option  to  select  “L  -  Log  Name”  to  change  the  name 
of  the  simulation  log,  “O  -  Output"  to  change  the  name  of  the  simulation  output  (data)  file, 
“R  -  Run  Simulation”  to  execute  the  backscatter  simulation,  or  “Q  •*  Quit”  to  return  to  the 
main  menu.  Type  the  appropriate  letter  and  press  ENTER. 

4.10  Changing  the  System  Configuration 

Selecting  “C  --  Configuration”  from  the  Main  Status  Menu  will  display  the  Master  Config¬ 
uration  Menu  as  shown  in  Figure  9.  Each  user  will  normally  set  up  a  configuration  tailored  to 
the  host  system,  applications,  and  configuration  preferences.  This  configuration  should  be  saved 
for  use  during  subsequent  sessions,  so  that  it  will  not  be  necessary  to  reset  the  configuration  from 
the  default  each  time  that  the  program  is  used.  The  features  that  make  up  a  given  configuration 
are  described  below. 

4.10.1  “N  -  Name”  Option 

This  option  prompts  the  user  for  a  new  configuration  name.  Type  in  a  name  of  eight 
characters  or  less  and  press  ENTER.  The  name  entered  the  user  immeuiaiely  becomes  the 
name  of  the  current  configuration. 
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N«m  of  Configuration  DEFAULT 
Propagation  Modal  ATMOSPHERIC  MODEL 
Raylolgh  Scattaring  Moda  ON 


N  —  Nana  L  —  Load  D  —  Display 

R  --  Run  Sat  P  —  Propagation  Moda  S  —  Scattaring  Moda 

0  ••  Quit 

I  Entar  Selection,  or  Praas  ENTER  Kay  for  Default  (Highlighted)  Selection. 

i _ 

Figure  9.  Muter  Configuration  Menu 

4.10.2  “L  -  Load”  Option 

This  option  converts  the  configuration  name  to  a  file  name.  If  the  file  exists,  the  program 
then  loads  the  configuration  information  contained  in  the  file.  If  the  file  does  not  exist,  the 
current  configuration  information  is  retained  and  referenced  by  the  new  name. 

4.10.3  “D  -  Display”  Option 

This  option  displays  the  Terminal  Display  Submenu,  which  is  shown  in  Figure  10.  The 
terminal  display  parameters  that  are  set  from  this  submenu  enable  the  program  to  adapt  the 
menu  control  system  to  the  characteristics  of  the  user’s  terminal. 

Changes  to  the  display  characteristics  that  are  selected  from  this  submenu  will  take  effect 
upon  quitting  the  Muter  Configuration  Menu  unless  the  ”C  -  Cancel”  option  is  selected  from 
the  quit  submenu.  This  behavior  allows  the  user  to  correct  a  mistake  or  cancel  a  series  of 
inadvertent  changes  without  the  screen  becoming  jumbled. 
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T  --  Terminal 
SCREEN  «IZC: 


C  —  Control  Char  A  —  Atmos  Displ  Q  —  Quit 
L  —  Length  U  —  Width  M  —  Markina 


Cntar  Selection,  or  Preee  COTTA  Key  for  Dafault  (Hiflh lighted)  Selection. 


Figure  10.  Terminal  Display  Submenu 

4.10.3.1  “T  --  Terminal’’  Qpaga 

This  selection  displays  the  choices  for  terminal  type.  The  terminal  selection  consists  of  only 
two  options,  ANSI  or  GENERIC.  Users  of  IBM  or  IBM-compatible  computers  (such  as  the 
Zenith  248)  or  users  with  terminals  conforming  to  ANSI  standard  (VT-100,  etc.)  should  select 
the  "A  -  ANSI"  option  and  select  the  ANSI  control  character  sequences  from  the  Control 
Character  submenu  (see  Subsection  4.10.3.2.)  This  use  of  this  option  requires  the  use  of  the  PC- 
DOS  or  MS-DOS  ANSI.SYS  display  driver  which  must  be  called  from  the  AUTOEXEC.BAT 
file  during  computer  system  "boot-up."  For  systems  that  do  not  conform  to  the  ANSI  standard 
for  terminal  control,  the  "G  -  Generic"  option  must  be  selected. 

Selecting  the  ANSI  terminal  option  allows  the  program  to  use  ANSI  terminal  control  ("es¬ 
cape")  strings  for  screen  control.  When  the  ANSI  option  is  in  effect,  the  following  features  are 
used  in  the  menu  displays: 

•  Scrolling  is  eliminated  by  removing  irrelevant  information  from  the  screen  and  positioning 

new  information  where  it  is  desired. 
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•  Default  menu  selections  are  highlighted  by  reverse  video.  (The  generic  display  driver  writes 

a  separate  line  to  indicate  the  default  selection.) 

•  Program  messages  signifying  normal  execution  are  displayed  in  high  intensity  text. 

-  Error  messages  are  displayed  in  blinking  high  intensity  text. 

These  enhancements  provide  a  more  intuitive  and  less  distracting  interface  for  most  users,  but 
may  not  be  available  on  all  hardware.  The  terminal  must  support  the  ANSI  protocol,  and  the 
computer  must  be  capable  of  transmitting  two  control  characters  described  below  to  the  terminal. 
All  of  the  menus  and  submenus  displayed  in  this  users  guide  were  generated  on  an  ANSI  terminal. 

4.10.3.2  “C  -  Control  Char* *  Potion 

The  Control  Character  Submenu,  which  is  shown  in  Figure  11,  is  displayed  by  selecting 
the  “ C  *  Control  Char”  option  from  the  Terminal  Display  submenu.  This  submenu  is  used  to 
define  the  parameters  used  to  determine  the  “escape”  control  character  and  the  “set  graphics 
rendition”  character,  in  BACKSCAT  parlance,  which  is  used  to  activate  and  deactivate  reverse 
video,  high  intensity,  blinking,  and  other  display  attributes. 

Since  many  computer  systems  support  either  of  two  standard  character  sets  (the  8-bit  ANSI 
set  and  IBM’s  EBCDIC  set),  the  control  characters  from  these  character  sets  are  built  into  the 
program  and  can  be  set  by  selecting  the  “A  -  ANSI”  and  “I  -  IBM”  options,  respectively. 
Fr*r  system  configurations  not  conforming  to  the  either  of  these  standards,  the  “E  -  ESCAPE” 
and  ‘‘vj  -  Graphics”  options  may  be  used  to  set  strings  to  be  substituted  for  the  escape  and  set 
graphics  rendition  characters  respectively.  Up  to  five  characters  may  be  entered  by  their  decimal 
positions  in  the  collating  sequence  of  the  host  computer  for  each  of  these  strings.  This  entry 
format  makes  the  entire  character  set  of  the  host  computer  available. 

410.3.3  VA-  AunoLPigpr.’  Option 

This  option  displays  the  choices  regarding  how  the  information  about  the  atmospheric 
aerosols  will  be  displayed  on  the  screen.  If  “N  ••  Normal”  is  selected,  the  information  about 
the  highest  aerosol  layers  in  the  atmosphere  will  appear  near  the  top  of  the  screen  and  that  about 
the  lowest  near  the  bottom  (see  Figure  7.)  If  “I  -  Inverted”  Is  selected,  this  sequence  will  be 
reversed.  Type  N  or  I  and  then  press  ENTER. 
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ESCAPE  Saquanca  }7 

ATTRIS  Saquanca  709 


i 

L 

Figure  11. 

4.10.3.4  Screen  Sire  Qpripni 

The  “L  -  Length"  and  “W  -  Width"  option*  respectively  let  the  number  of  line*  (verti¬ 
cally)  and  number  of  characters  (horizontally)  of  the  screen.  The  "M  -  Margins"  option  sets  the 
width  of  the  left  and  right  margins,  in  characters,  in  the  status  display  area.  Type  the  appropriate 
response  and  press  ENTER.  On  systems  that  reserve  one  or  more  lines  at  the  top  or  bottom  of 
the  screen  for  sums  display,  the  screen  length  should  be  set  to  the  number  of  availablt  lines. 

4.10.3.5  "0  -  Quit"  Potion 

This  selection  returns  control  to  the  Master  Configuration  Menu.  Type  Q  then  press  ENTER. 


ANSI  CONTROL)  C  --  Eacapa  Saquanca 

A  —  ASCII  standard  Sat  I  —  XSK  (MCDIC)  Sat 


e  —  crapnica  saquanca 
Q  --  Quit 


Entar  Salaction,  or  Rraaa  ENTER  Kay  tor  Datauit  (Hiqhllqhtad)  Salaotion. 


The  Control  Character  Submenu 
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4.10.4  "R  -  Run  Set”  Option 

This  option  displays  the  current  run  set,  as  shown  in  the  Run  Set  Submenu  in  Figure  12, 
consisting  of  die  lidar  system  name,  the  atmosphere  name,  the  propagation  profile  name,  the 
molecular  absorption  profile  name,  the  simulation  log  name,  and  the  simulation  output  name. 
The  nin  set  submenu  allows  each  of  these  to  be  changed.  This  submenu  is  included  primarily 
as  a  convenience  for  changing  the  default  run  set  elements  to  be  saved  in  a  configuration  file. 
Each  element  of  the  run  set  can  also  be  selected  from  another  menu  outside  the  control  of  the 
Master  Configuration  Menu. 


LI  DAK  Systea 

DEfAULT 

Atmospheric  Model 

DEFAULT 

Propagation  Profile 

TEST 

Molecular  Absorption 

(HONE) 

Simulation  Log 

DEFAULT 

Siaulatlon  Output 

DEFAULT 

PUN  8tTi 

—  Atmospheric  Model 

—  Molecular  Absorption 
--  Siaulatlon  Output 

S  --  LIDAR  Systea 
r  --  Propagation  Profile 
L  —  Siaulatlon  Log 

Q  --  Cult 

i  enter  selection,  or  frees  EMTER  Key  for  Default  (Highlighted)  Selection. 


Figure  12.  The  Run  Set  Submenu 


440.5  “P  -  propagation  Mode”  Option 

This  selection  displays  the  options  that  allow  the  user  to  select  either  of  two  modes  by  which 
the  backscauer  can  be  computed.  If  .the  “A  —  Atmospheric  Model”  option  is  selected,  the 
simulation  will  compute  the  atmospheric  propagation  coefficients  from  the  AFGL  atmospheric 
model  and  write  them  to  the  file  associated  with  the  propagation  profile  name.  If  the  “P  - 
Propagation  Profile”  option  is  selected,  the  simulation  will  read  the  propagation  coefficients 
from  the  file  associated  with  the  propagation  profile  name.  (2>ectipn  4.6  describes  the  format 
requirements  of  a  user-supplied  propagation  profile.) 

4.10.6  “S  -  Scattering  Mode”  Option 

This  selection  displays  the  options  regarding  the  use  of  Rfylejgfe  featuring.  flip  #ser  fp*y 
select  ‘*y  *•  yes”  to  include  Jfciyleigb  scattering  or  "N  -  No”  tp  neglect  Rayleigh  scattering. 
The  status  of  the  scattering  mode  is  indicated  as  “ON"  if  Raylpigh  scattering  is  included  or 
“OJFF”  if  it  is  not. 

4.10.7  "Q  -  Quit”  Option 

This  option  exits  the  Master  Configuration  Menu  and  returns  control  tp  the  Kfaip  Status 
Menu.  To  select  this  option,  type  Q  and  press  ENTER. 

4.11  Execution  in  Batch  Mode 

The  BACKSCAT  program  may  be  executed  in  batch  mode  by  redirecting  standard  input 
to  a  script  file  using  the  redirection  capabilities  of  the  operating  system.  (Under  PC-DOS 
and  MS-DOS,  the  <  and  >  characters  are  used  to  redirect  standard  input  and  standard  output 
respectively.)  A  script  file  is  a  input  text  file  containing  the  consecutive  entries  that  would  be 
supplied  interactively  through  the  menu  system;  it  can  be  prepaid  using  an  ASCII  te$t  fldHpr. 
All  functions  supported  by  the  menu  system  can  be  performed  through  appropriate  entries  jn  the 
script  file.  The  script  file  shown  in  Figure  13  loads  a  configuration  renames  die  ptjfput,  executes 
the  simulation,  and  exits  the  program;  the  parenthetical  continent*  are  optional.  All  fegt  jn  a 
script  file  is  automatically  converted  to  UPPER  CASE  in  the  ifme  manner  a#  interactive  input. 

For  batch  mode  execution,  standard  output  should  also  be  redirected  to  a  “null”  or  "phan¬ 
tom”  file  or  device  that  serves  to  suppress  output  (for  example,  the  PC-DOS  or  Mg-QQS  device 
“NUL”)  if  possible,  as  standard  output  is  used  only  to  display  the  screens  and  messages  of  the 
menu  control  system.  In  case  of  difficulty,  this  information  can  instead  be  redirected  to  a  file 
and  subsequently  matched  with  the  script  file  as  a  diagnostic  measure. 


24 


O  (A 


KY COIF 10 


HYOUTPOT 

R 

Q 

Q 


(initial  configuration  name) 

(select  simulation  from  main  menu) 

(select  output  name  from  simulation  menu) 
(new  output  name) 

(select  run  simulation) 

(quit  simulation  menu) 

(quit  program) 


Figure  13.  A  Sample  Script  File  for  Batch  Mode  Execution  of  BACKSCAT 


5  SAMPLE  OUTPUT 


The  simulation,  whether  executed  interactively  or  in  batch  mode,  always  generates  an  output 
file  and  a  log  file.  The  lidar  output  file  contains  the  backscatter  data  in  standard  tabular  form  for 
input  into  any  standard  plotting  program;  many  such  programs  are  available  commercially.  The 
columns  of  the  output  file  are:  (1)  range  (kilometers),  (2)  altitude  (kilometers),  (3)  cumulative 
attenuation  factor  (“optical  depth"),  (4)  lidar  backscatter  (Watts),  (3)  normalized  backscatter, 
(6)  range-compensated  backscatter  (Watt-square  meters),  and  (7)  normalized  range-compensated 
backscatter  (square  meters). 

The  log  file  provides  a  trace  of  the  simulation  execution  that  can  be  used  to  validate  results. 
Log  files  are  self-documenting,  and  their  contents  varies  depending  upon  the  simulation  mode 
options  that  are  in  effect  for  the  particular  execution.  A  sample  log  file  is  shown  In  Figure  14 
with  only  a  partial  prim  out  of  the  lidar  backscatter  return  data. 

When  the  propagation  mode  is  selected  to  the  atmospheric  model,  the  simulation  also  gen¬ 
erates  a  propagation  profile  file.  The  propagation  profile  file  is  also  in  standard  tabular  form  for 
input  into  any  standard  plotting  program;  it  can  additionally  be  input  into  he  simulation  on  sub¬ 
sequent  executions  by  selecting  the  propagation  mode  to  the  propagation  profile.  The  columns 
of  the  propagation  profile  are:  (1)  altitude  (kilometers),  (2)  aerosol  extinction  coefficient  (in¬ 
verse  kilometers),  (3)  aerosol  scattering  coefficient  (inverse  kilometers),  (4)  aerosol  absorption 
coefficient  (inverse  kilometers),  (3)  aerosol  backscatter  coefficient  (inverse  meten-steradians), 
(6)  molecular  extinction  coefficient  (inverse  kilcmeten),  and  (7)  molecular  backscatter  coeffi¬ 
cient  (inverse  meters-steradians). 

Figures  13-17  show  the  results  from  a  sample  run.  Figure  13  shows  the  propagation  profile 
created  by  the  specific  choice  of  atmospheric  properties  and  aerosols.  The  scattering  coefficient 
is  hidden  behind  the  total  extinction  curve  since  the  absorption  at  this  wavelength  is  relatively 
small.  The  backscatter  coefficient  is  also  shown  on  the  same  graph  with  the  same  scale  but 
different  units.  Figure  16  and  17  show  the  lidar  backscatter  return  with  the  1/J?*  dependence 
included  and  removed,  respectively. 
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LIDAR  Backacatter  Simulation 
(Version  1.1) 

Simulation  Log 


LIDAR  SYSTEM 


Transmitter  Parameters: 

Wavelength  (urn) 
Pulse  Energy  (J) 
Pulse  Duration  (us) 


Receiver  Parameters: 

Aperture  Diameter  (cm) 
Obscuration  Diameter  (cm) 


Signal  Processor  Parameters: 

Nearest  (Minimum)  Range  (km)  0.0 

farthest  (Maximum)  Range  (km)  100.0 

Range  Resolution  (Kilometers)  .5 


Viewing  Aspect: 

Sensor  Height  (km) 

viewing  Azimuth  Angle  (Deg) 

Viewing  Elevation  Angle  (Deg) 


ATMOSPHERIC  MODEL 


Boundary  Layer: 

Aerosol 

Visibility  (km) 

Relative  Humidity  (Percent) 


RURAL 

23 

70.0 


Troposphere: 

Relative  Humidity  (Percent) 


70.0 


Stratosphere: 

Aerosol  Extinction 
Aerosol  Profile 


BACKGROUND 

BACKGROUND 


Upper  Atmosphere: 

Aerosol  Level 


NORMAL 


Aerosol  Distribution: 


FALL/ WINTER 


Layer  Heights: 

Boundary  Layer  (km) 
Troposphere  (km) 
Stratosphere  (km) 

Top  of  Atmosphere  (km) 


2.0 

9.0 

29.0 

100.0 


Propagation  profile  written  to  file  TEST.PFL 

Figure  14.  Sample  Log  File  Output 
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PROPAGATION  PROFILt 


Height 

Satinet 

Scatter 

Abaorp 

B'acattar 

Hoi  Ixt 

Mol  B'acat 

Out) 

U/Aft) 

(1/kn) 

(1/hRv) 

(1/n-ar) 

(l/ka) 

(1/at-ar) 

.00 

1 . 5601-01 

1 . 494E-01 

8.491E-03 

3 . 466E-06 

1.1601-02 

1.3851-06 

1.00 

9.910E-02 

9.373E-02 

5 . 3251-03 

2 . 174E-06 

1.0601-02 

1.265t-06 

1.50 

7.920E-C? 

7 . 491E-02 

4 . 256E-03 

1.737E-06 

1.0001-02 

1.1941-06 

2.00 

6.2101-02 

5.874E-02 

3.337E-03 

1 . 362E-06 

9.5501-03 

1.1401-06 

2.00 

2.7201-02 

2 . 621C-02 

9 . 931E-04 

5 . 918E-07 

0.6201-03 

1.0291-06 

4.00 

1.2001-02 

1 . 1S6E-02 

4 . 381E-04 

2 . 611E-07 

7.7701-03 

9.2751-07 

$.00 

4.6501-03 

4.673E-03 

1.771E-04 

1 . 055E-07 

C.  9901-03 

6.3441-07 

«.00 

3.540C-03 

3.411E-03 

1.292E-04 

7 . 702E-08 

6.2601-03 

7.4721-07 

7 . 00 

2 . 300E-03 

2.216E-03 

6 . 397E-0S 

5 . 004E-06 

5.6001-03 

6.684E-07 

e.oo 

1.4101-03 

1.359E-03 

S.146E-05 

3 . 068E-08 

4.9901-03 

5.9561-07 

9.00 

9.600E-04 

9 . 442E-04 

3.576E-05 

2 . 132E-08 

4.4301-03 

5.2661-07 

10.00 

7.670E-04 

7.870E-04 

4 . 658E-1 1 

1.124E-08 

3.9201-03 

4.6791-07 

11.00 

7 . 140E-04 

7.140E-04 

4 . 226E-11 

1 . 020E-06 

3.460t-03 

4.1301-07 

12.00 

6 . 630E-04 

6.630E-04 

3 . 924E-11 

9 . 468E-09 

2.9601-03 

3.5331-07 

13.00 

6.220E-04 

6.220E-04 

3 . 682E-11 

8 . 662E-09 

2.5301-03 

3.0201-07 

14.00 

6 . 450E-04 

6 . 450E-04 

3.618E-11 

9 . 211E-09 

2.1601-03 

2.5781-07 

15.00 

6 . 430E-04 

6 . 430E-04 

3. 606E-11 

9 . 162E-09 

1.8501-03 

2.2081-07 

16.00 

6 . 410E-04 

6 . 410E-04 

3.794E-11 

9. 153E-09 

1.5601-03 

1.8661-07 

17.00 

6 . 010E-04 

6 . 010E-04 

3.5S7E-11 

8 . 562E-09 

1.3501-03 

1.6111-07 

16.00 

5.630E-04 

5 . 630E-04 

3.332E-11 

6.040E-09 

1.1501-03 

1.3731-07 

19.00 

4 . 920E-04 

4 . 920E-04 

2 . 912E-11 

7 . 026E-09 

9.6701-04 

1.178t-07 

20.00 

4.230E-04 

4.230E-04 

2.504E-11 

6 . 040E-09 

6.4401-04 

1.0071-07 

21.00 

3.520C-04 

3 . 520E-04 

2 . 064E-11 

5.027E-09 

7.1101-04 

6.5701-00 

22.00 

2.960E-04 

2 . 960E-04 

1 . 752E-11 

4.227E-09 

6.1201-04 

7.3051-08 

23.00 

2 . 420E-04 

2 . 420E-04 

1 . 432E-11 

3. 4561-09 

5.2201-04 

6.2311-06 

24.00 

1 . 9001-04 

1 . 900E-04 

1. 125E-11 

2.713E-09 

4.4501-04 

5.3121-08 

25.00 

1 . 500E-04 

1 . 500E-04 

S.879E-12 

2.1421-09 

3.6001-04 

4.5361-06 

26.00 

1 . 1S0E-04 

1 . 150E-04 

6 . 607E-12 

1.6421-09 

3.2501-04 

3.8791-06 

27.00 

6 . 950E-05 

6.950E-05 

5 . 298E-12 

1.2761-09 

2.7601-04 

3.3161-06 

26.00 

6.700E-05 

6.700E-05 

3.966E-12 

9.566E-10 

2.3601-04 

2.8411-08 

29.00 

5.200E-05 

5.200E-0S 

3.078E-12 

7 . 426E-10 

2.0401-04 

2.4351-06 

30.00 

3 . 320E-05 

3.320E-05 

1 . 965E-12 

4.741E-10 

1.7501-04 

2.0891-06 

35.00 

1 . 650E-05 

1 . 650E-05 

9.767E-13 

2. 3561-10 

6.0301-05 

9.5651-09 

40.00 

6.000E-06 

6 . 000E-06 

4.735E-13 

1.142E-10 

3.7901-05 

4.5241-09 

45.00 

4.020E-06 

4 . 020E-06 

2.379E-13 

5.7411-11 

1.6701-05 

2.2321-09 

50.00 

2 . 100E-06 

2 . 100E-06 

1 .243E-13 

2.9991-11 

9.7401-06 

1.1631-09 

55.00 

1 . 090E-06 

1 . 090E-06 

6 . 452E-14 

1.SS7E-11 

5.2001-06 

6.2071-10 

60.00 

5.760E-07 

5 . 760E-07 

3 . 421E-14 

6.2541-12 

2.7201-06 

3.2471-10 

65.00 

3.050E-07 

3 . 050E-07 

1 . 605E-14 

4.3551-12 

1.4801-06 

1.7671-10 

70.00 

1 . 600E-07 

1 . 600E-07 

9.471E-15 

2.265E-12 

6.0501-07 

9.6091-11 

75.00 

6 . 9S0E-06 

6 . 950E-06 

4 . 114E-15 

9.9251-13 

3.3001-07 

3.9391-11 

60.00 

2 . 900E-08 

2.900E-08 

1.717E-15 

4.1411-13 

1.3501-07 

1.6111-11 

65.00 

1.200E-06 

1.200E-08 

7 . 103E-16 

1 .7141-13 

5.6501-06 

6,7441-12 

90.00 

5.100E-09 

5.100E-09 

3 . 019E-16 

7. 2831-14 

2.4001-06 

2.6651-12 

95.00 

2 . 150E-09 

2 . 150E-09 

1.273E-16 

3.0701-14 

1.0001-06 

1.1941-12 

100.00 

9.300E-10 

9 . 300E-10 

5.505E-17 

1.3261-14 

4.2001-09 

5.0131-13 

Molecular  acattaring  mod*!  la  ON. 

No  abaorption  from  molacular  raaonaneat. 
Simulation  output  written  to  file  TK9T.OAT 


Figure  14.  Sample  Log  File  Output  (cont’d) 
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LXDAR  EACKSCATTER 


Ranqa 

(Km) 

5.000E-01 
l.OOOE+OO 
1 . SOOE+OO 

2 .  OOOE+OO 
2 . JOOE+OO 
3.000E+00 
3 . 500E+-00 
4  .  OOOE+OO 
4 . 500E+00 
5 . OOOE+OO 

5. $00E*00 

6 .  OOOE+OO 
€ . SOCE+OO 
7.000E+00 

7 .  500E+00 

e.oooE+oo 

8 . SOOE+OO 
9 . 000E+00 
9 . 50CE+00 

1 .  000E+01 

1  .OSCE+Ol 
1 . 100E+01 
1 . 150E4-01 
1.200E+01 
1 . 250E+01 
1 .300E+01 
X . 350E+01 
1 . 400E+01 
1.4S0Et>0l 
1.500E+01 

1 . SSQE+Ol 
1 . 600E+01 
X . 680E+01 
X . 700E+01 
1.7SOE+01 
1 . 800E+01 
1.850E+C1 
1.900E+0X 
X . 9S0E+01 

2 .  OOOE+Oi 

2 . 050E+01 
2 .  100E+01 
2 . XSOE+OX 
2 . 200E+01 
2.280E+01 
2 , 300E+01 
2 . 350E+01 
2.400E+01 
2 .  4S0E+O1 
2 . 500E+01 


Haight 

(Km) 

3 .  S36E-01 
7.071E-01 
x  .oeiE+oo 
X . 414E+00 
X.789E+00 
2 . 121E+00 
2.475E+00 
2 . 628E+00 
3.X82E+00 
3 . 536E+00 

3 . 889E+00 
4 . 20E+00 
4 . 596E+00 
4 . 950E+00 
5 . 303E+O0 
5 . 857E+00 
6 . OlOE+OO 
6.364E+00 
6.718E+0C 
7 . 071E*00 

7.425E*00 
7.778E+00 
8 . 132E+00 
8.48SE+00 
8.839EfOO 
9 .  192E+00 
9.S46E+00 
9 . 899E+00 
X.028E+01 
1 . 061E+01 

1.098E+01 
X .  X3XE+01 
1.X67E+01 
X  .202E+01 
1 .237E+01 
1.273E+01 
1.308E+01 
X .  344E+01 
1 . 379E+01 
1.414E+01 

X .  4S0E+01 
X .  48SE+01 
1 .520E+01 
X . 5S8E+0X 
X . S91E+01 
X . 826E+01 
X . 662E+01 
1.697E+01 
X  .732E+01 
1 .768E+0X 


Optical 

Dapth 

(-) 

7 . 9S1E*02 
1 . 484E-0X 
2.070E-01 
2.570E-01 
3.001E-0X 
3.369E-01 
3.874E-01 
3.916E-01 
4.103E-01 
4.253E-01 

4.375E-01 
4.473E-01 
4.556E-01 
4.824E-01 
4 . 683E-0X 
4 .737E-01 
4.788E-01 
4.835E-C1 
4.879E-01 
4 . 919E-01 

4.957E-01 
4 . 992E-01 
5 . 024E-01 
S.0J5E-01 
S.084E-01 
5.1HE-01 
S.136E-01 
5.181E-01 
5.X84E-01 
5.207E-01 

5.228E-01 
5 . 249E-01 
5.268E-01 
S.287E-0X 
S.304E-01 
5.321E-01 
S.337E-01 
5 . 352E-01 
5 . 367E-01 
5.381E-01 

5.395E-OX 

5 .  408E-OX 

5.420E-01 

5.432E-01 

5 . 444E-01 

S.4S5E-01 

5 . 483E-01 

5.475E-01 

5.485E-01 

5 . 494E-01 
• 


Udar 

Raturn 

(«) 

1 .748E-03 
3.370E-04 
X . 168E-04 
5 . 3X2E-0S 
2 .  7  92E-05 
1.597E-05 
9.603E-06 
5.958E-06 
3.946E-06 
2.794E-06 

2 . 006E-06 
X . 508E-06 
X . X  65E-06 
9.075E-07 
7.423E-07 
i . 1 62E-07 
5 . 148E-07 
4 .342E-07 
3.679E-07 
3 . 134E-07 

2.692E-07 
2.319E-07 
2 . 008E-07 
1 .753E-07 
1.533E-07 
1.345E-07 
1.182E-07 
1.040E-07 
9.232E-08 
8 . 223E-08 

7 . 327E-08 
6.507E-08 
S.773E-08 
5 .  118E-08 

4 .  S69E-08 
4.074E-08 
3 . 635E-08 
3 . 263E-08 
2 . 924E-08 
2 . 626E-08 

2.368E-08 
2 .  X32E-08 
X .  92SE-06 
X.740E-08 
1.571E-08 
1.423E-08 
X.290E-08 
X .  X66E-06 
X .  060E-06 
9.816E-09 


Nomalitad 

Raturn 

(-) 

1.748E-09 

3.370E-10 

1 .  168E-10 
5.312E-XX 

2 .  792E-11 
1.597E-1X 
9 . 603E-X2 
5.958E-12 
3 . 94CE-X2 
2.794E-12 

2 . 006E-12 
1 . 508E-12 
1.165E-12 
9.075E-13 
7.423E-13 
6 . 162E-13 
S.148E-13 
4 . 342E-13 
3.679E-13 
3.134E-13 

2.692E-13 
2 . 3X9E-X3 
2 . 008E-13 
X.7S3E-13 
X.533E-13 
1.345E-X3 
X .  182E-X3 
X .  040E-13 
9.232E-X4 
8 . 223E-14 

7.327E-14 
6.507E-14 
5.773E-14 
5.118E-14 
4 . 569E-14 
4.074E-14 
3 . 635E-14 
3.263E-14 
2.924E-X4 
2.626E-X4 

2.388E-X4 
2.132E-14 
X .  925E-14 
X.740E-X4 
1 . 571E-14 
X . 423E-14 
1 . 290E-14 
X . 166E-14 
X . 060E-14 
9.816E-15 


Rng.  Xnd. 

Lidar 
(W-aq  m) 

4 . 369E+02 
3 . 370E+02 
2 . 628E+02 
2 . 125E+02 
1.743E+02 
1 . 437E+02 
X .  178E+02 
9 . 533E+01 
7 . 991E+01 
6 . 986E+01 

S . 068E+01 
5 . 429E+01 
4 . 924E+01 
4 . 447E+01 
4.176E+01 
3 . 944E+01 
3.720E+01 
3 . 5X7E+01 
3 . 321E+01 
3 . 134E+01 

2.968E+01 
2 . 8C6E+01 
2 . 6S6E+01 
2.52SE401 
2 . 395E+01 
2.273E+01 
2 . 15SE+01 
2 . 039E+01 
X  .  941E+01 
X . 850E+01 

1.760E+01 
X  .  666E+01 
1.572E+0X 
1.479E+01 
X .  399E+0X 
X . 320E+01 
1 . 244E+0X 
X.178E+01 
1.112E+01 
X  .  050E+01 

9.953E+00 
9 . 404E+00 
8 . 898E+00 
8.424E+00 
7 . 952E+00 
7 , 528E+00 
7 . 121E+00 
8.717E+00 
( ,  360E+00 

6 . 010E+00 
• 


Figure  14.  Sample  Log  File  Output  (corn’d) 


Normallted 
Rat urn 
<aq  m) 

4 . 369E-04 
3.370E«04 
2 . 628E-04 
2.125E-04 
X.745E-04 
1 . 437E~04 
1.176E-04 
9 .533E-05 
7 . 991E-0S 
8 . 986E-05 

6.068E-05 
5 . 429E-0S 
4.924E-05 
4.447E-0S 
4.176E-US 
3 . 944E-05 
3.720E-0* 
3.517E-05 
3 . 321E-05 
3 . 134E-05 

2 . 968E-0S 
2 . 806E-05 
2 . 858E-0S 
2 . 525E-05 
2 . 395E-05 
2.273E-05 
2 . 155E-05 
2 . 039E-05 
X .  941E-05 
X . 850E-0S 

X.780E-05 
1.666E-0S 
1 . 5"2E-0S 
1.479E-0S 
1 . 399E-05 
X . 320E-05 
1.244E-05 
1.178E-05 
1 . 112E-05 
1.050E-0S 

9.953E-06 
9 . 404E-06 
8.898E-06 
8 . 42  4E-06 
7 . 9S2E-06 
7.528E-08 
7 . 121E-06 
8.717E-06 
6.360E-06 
8.0X0E-04 
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Attenuation  Coefficients  (1/km) 


Altitude  (km) 


Figure  15.  Atmospheric  Attenuation  and  Backscatter  Coefficients  Versus  Altitude 
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Normalized  Return 


Normalized  Return  (w/o  R2) 


10 

10 

10 

10 

10 


0.0  40.0  80.0  120.0 

Altitude  /  Range  (km) 


Figure  17.  Range  Impendent  Lidar  Backscatter  Return 


6  SUMMARY  AND  RECOMMENDATIONS  FOR  FUTURE  MODIFICATIONS 

BACKS  CAT  simulates  the  return  form  a  backscatter  lidar  system.  The  atmospheric  and 
lidar  system  conditions  that  define  the  nature  of  the  return  are  under  operator  control. 

This  version  of  BACKSCAT  represents  a  first  attempt  to  simulate  a  complete  lidar  system. 
A  number  of  improvements  to  BACKSCAT  could  be  made  to  make  it  more  representative  of 
an  actual,  fieldable  system. 

6.1  Improvement  of  Model  Physics 

6.1.1  Multiple  Pulses 

BACKSCAT  simulates  the  return  from  a  single  laser  pulse.  BACKSCAT  could  be  modified 
to  include  multiple  pulses.  Different  pulse  averaging  techniques  could  also  be  included. 

6.1.2  Addition  of  Noise 

There  are  irregularities  between  lidar  pulses  that  are  not  accounted  for  in  BACKSCAT. 
These  differences  can  be  accounted  for  through  the  addition  of  simulated  noise.  Depending  on 
the  degree  of  sophistication  desired,  one  could  include  one  generic  noise  term  or  separate  terms 
for  laser  and  receiver  electronics  noise. 

6.1.3  Allowance  for  Atmospheric  Fluctuations 

BACKSCAT  assumes  a  static  atmosphere.  BACKSCAT  could  be  coupled  to  either  a 
database  of  time-dependent  atmospheric  values  or  to  an  external  program  that  generates  time 
varying  atmospheric  profiles. 

6.1.4  Consideration  of  Wind  Dependent  Aerosols 

The  aerosol  models  included  in  BACKSCAT  are  also  static  models.  The  wind  dependent 
aerosol  models  such  as  the  desert  and  Navy  maritime  models  should  be  added  to  the  list  of 
aerosols  that  can  be  considered. 

6.1.5  Consideration  of  Variable  Boundary  Heights 

The  boundary  heights  used  to  diferentiate  one  region  of  aerosols  from  another  are  currently 
fixed  In  BACKSCAT.  One  can  use  other  boundary  heights  but  that  would  entail  the  use  use  of  a 
customized  set  of  scaling  factors  (see  section  4.6).  It  is  recommended  that  the  code  be  modified 
to  allow  the  use  of  variable  boundary  heights.  This  modification  would  entail  the  addition  of 
extrapolation  routines  that  would  automatically  adjust  the  scaling  factors  with  the  new  boundary 
heights. 
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6.1.6  Alternate  Calculation  of  Rayleigh  Scattering 

BACKSCAT  currently  use  a  default  profile  for  Rayleigh  scattering.  The  code  should  be 
modified  to  allow  the  calculation  of  a  Rayleigh  scattering  profile  via  the  use  of  radiosonde  data 
or  alternate  model  atmospheres. 

6.1.7  Variable  Molecular  Absorption  Profiles 

BACKSCAT  includes  a  default  molecular  absorption  profile  (see  section  4.8)  based  on  a 
midlatitude  atmosphere.  BACKSCAT  should  be  modified  to  allow  the  use  of  other  molecular 
absorption  profiles,  such  as  those  that  can  be  produced  by  FASCODE.  In  a  further  modification, 
BACKSCAT  could  be  coupled  to  a  special  version  of  FASCODE  to  calculate  the  molecular 
absorption  profile  for  given  laser  lines. 

6.1.8  Allow  for  Scattering  at  Other  Than  180  Degrees 

BACKSCAT  is  presently  set  up  to  model  backscattering  at  180  degrees.  The  code  should 
be  modified  to  allow  for  experimental  configurations  in  which  the  transmitter  and  receiver  are 
at  different  locations.  This  modification  would  require  the  use  of  a  larger  data  base  of  phase 
functions  and  separate  atmospheric  path  definitions  for  the  transmitter  and  receiver  paths.  This 
modification  would  also  require  the  addition  of  field-of-view  effects  in  the  processor  module. 

6. 1 .9  Addition  of  Field-of-View  Effects 

Field-of-view  effects  are  not  now  currently  included.  In  order  to  adequately  simulate  a 
fielded  system,  they  should  be  incorporated.  As  noted  in  the  above  recommendation,  these 
effects  would  have  to  be  included  if  scattering  at  other  than  180  degrees  was  included. 

6.2  User  Features 

6.2.1  On-Screen  Help  Function 

BACKSCAT  does  not  include  any  built-in  help  capability  in  which  the  user  can  press  a  key 
and  get  information  about  particular  features.  The  addition  of  an  on-screen  help  system  would 
aid  the  user  in  the  use  of  BACKSCAT.  This  feature  would  be  intended  to  augment  the  user’s 
manual,  not  replace  it 

6.2.2  Alternate  Units 

BACKSCAT  assumes  a  fixed  set  of  units  for  the  various  parameters.  The  user  should  be 
given  the  option  of  using  other  units  for  parameters,  such  as  the  laser  power. 
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6J  Computer  Code  Changes 

The  computer  code  was  written  in  Fortran  77  which  is  a  highly  suitable  language  for  scientific 
calculations  but  a  restrictive  language  for  developing  menu-based,  user  interface  input  systems. 
The  C  computer  language  is  more  suitable  for  the  development  of  menu-based  input  systems. 

It  is  recommended  that  the  user-interface  system  of  BACKSCAT  be  rewritten  in  C.  The 
portion  of  the  code  devoted  to  the  scientific  calculations  should  be  kept  in  Fortran  77.  This 
portion  of  the  code  would  be  connected  to  the  C- based  user  interface  udlizing  the  mixed-language 
capabilities  of  C. 
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APPENDIX  A  IM  of  Routines 

The  following  routines  comprise  she  BACKSCAT  backsotter  li4w  simulation  program.  The 
listing  describes  the  type  of  routine  and  a  brief  description  of  what  the  routine  does.  All  of  the 
routines  were  written  and  compiled  using  Microsoft  ^  FORTRAN  version  4.01  of  FORTRAN 
*77. 


Name _ Type  of  Routine _  P*"****0" 


BACKSCAT 

(Mein  Program) 

atmpfl 

LOGICAL  FUNCTION 

ATM  ST  A 

SUBROUTINE 

EDATMS 

SUBROUTINE 

EDCNfG 

SUBROUTINE 

EDUDR 

SUBROUTINE 

EXTINT 

LOGICAL  FUNCTION 

FBSCAT 

LOGICAL  FUNCTION 

GETATM 

SUBROUTINE 

GETCFG 

SUBROUTINE 

GETLDR 

SUBROUTINE 

GETMOD 

LOOICAL  FUNCTION 

GETPFL 

SUBROUTINE 

GETRES 

SUBROUTINE 

GETSCL 

LOOICAL  FUNCTION 

INTEGR 

SUBROUTINE 

INTRPL 

SUBROUTINE 

JJF1LE 

CHARACTER  •(•)  FUNCTION 

LDRCLT 

SUBROUTINE 

LIDAR 

SUBROUTINE 

MES2AG 

SUBROUTINE 

NEWNAM 

SUBROUTINE 

NEWSEQ 

SUBROUTINE 

PHGET 

LOOICAL  FUNCTION 

PROMPT 

SUBROUTINE 

PRSATM 

LOOICAL  FUNCTION 

READUC 

SUBROUTINE 

REALNO 

SUBROUTINE 

RHTNT 

SUBROUTINE 

Driver  for  main  menu  and  simulation  menu 
generates  propagation  profile  from  atmospheric  model 
Generates  master  atmospheric  status  display 
Edits  atmospheric  structure,  driving  master  atmospheric  struc¬ 
ture  menu  and  all  subordinate  menus 
Edits  program  configuration,  driving  master  configuration  menu 
and  all  subordinate  menus 

Frii>«  lidar  system,  driving  master  Udar  system  menu  and  all 
subordinate  menus 

Determines  extinction  coefficients  #  Udar  system  wavelength 
Calculates  backscaiter  normalization  factor  based  on  extinction 
and  visibility 

atmospheric  structure  from  file 
Loads  configuration  from  file 
Loads  lidar  system  from  file 
Reads  data  from  atmospheric  model  file 
Loads  propagation  profile  data 
Loads  molecular  resonance  absorption  data 
Reads  vertical  profile  scaling  data 
Prompts  user  for  integer  input  and  gets  value 
Determines  interpolation  point  in  index  array 
Converts  internal  name  to  file  name  for  specified  file  type 
Computes  lidar  attenuation  and  backacatter 
Converts  Udar  system  parameters  from  menu  interface  format 
to  simulation  format 

Displays  text  on  screen  message  line  with  specified  attributes 
Displays  prompt  and  gets  new  name 
Oets  new  ANSI  control  character  sequence 
Locates  correct  atmospheric  model  phase  function  and  Interpo¬ 
lates  on  relative  humidity  and  wavelength 
Displays  text  on  prompt  line  with  specified  attributes 
Panes  the  atmospheric  structure  specifications 
Reads  input  string  and  convene  to  upper  caw 
Displays  prompt  for  user  real  number  input  and  gets  value 
Performs  appropriate  interpolations  and  computes  extinction  es 


i  function  of  relative  humidity 
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Nam 


SAVATM 

SAVCFG 

SAVLDR 

SAVPFL 

SCLMOL 

SELECT 

STATUS 

VISINT 

WAVINT 


Type  of  Routine  Description 


SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
LOGICAL  FUNCTION 


Saves  atmospheric  structure  to  (lie 
Saves  configuration  to  file 
Saves  lidar  system  to  file 
Saves  propagation  profile  to  file 

Scales  molecular  backseat  ter  and  extinction  coefficients  by  ver¬ 
tical  concentration  profile 
Displays  menu  and  gets  user  selection 
Displays  information  in  status  window 
Interpolates  visibility  scaling  factor 
Scales  propagation  parameter!  for  wavelength 
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APPENDIX  B  Operating  System  Dependencies 

Portability  was  a  major  criteria  in  the  design  of  the  BACKSCAT  program;  hence,  computer 
and  operating  system  dependencies  have  been  eliminated  where  practicable  and  otherwise  iso¬ 
lated  in  a  minimal  number  of  low-level  routines.  One  problem  which  continues  to  remain  is 
differences  in  file  naming  conventions  between  different  operating  systems.  The  ANSI  stan¬ 
dard  for  FORTRAN  ’77  states  that  “a  file  may  have  a  name...  The  set  of  allowable  names  is 
processor-determined  and  may  be  empty."* 

The  problem  created  by  non-standardization  of  file  naming  conventions  has  been  resolved 
in  the  BACKSCAT  code  through  the  use  of  a  function  named  JJFILE  to  convert  the  names  of 
various  data  sett  in  the  simulation  to  file  names  whenever  files  are  read  or  written.  The  function 
subroutine  named  JJFILE  that  is  supplied  with  the  simulation  was  written  for  PC-DOS  and  MS- 
DOS  operating  systems,  and  will  also  work  satisfactorily  under  UNIX,  but  should  be  reviewed 
and  revised  if  necessary  whenever  the  simulation  is  ported  to  any  other  operating  system. 

The  function  JJFILE  accepts  two  arguments,  both  of  which  are  of  type  CHARACTER,  and 
the  function  itself  is  also  of  type  CHARACTER.  The  fust  argument  is  the  name  associated  with 
the  data  set;  under  PC-DOS  or  MS-DOS  (and  also  under  OS/2  and  UNIX),  it  may  contain  a  drive 
and/or  path  specifier  if  portability  of  data  Hies  is  not  an  important  consideration.  The  second 
argument  is  a  rile  type  from  Table  B-l.  The  file  type  is  always  four  characters  in  length,  and  the 
first  character  is  always  a  period.  The  version  of  JJFILE  supplied  with  the  program  scans  the 
first  argument  looking  for  a  period  indicating  that  a  file  extension  is  already  present;  if  none  is 
found,  it  concatenates  the  rile  type  as  a  file  extension.  A  lidar  system  named  DEMO  would  be 
read  from  or  written  to  a  file  named  DEMO.LDR,  while  an  atmospheric  structure  named  TEST 
would  be  read  from  or  written  to  a  rile  named  TEST.ATM. 

The  rile  types  in  Table  B-l  were  selected  based  on  several  criteria.  First,  they  have  been 
chosen  to  be  mnemonic;  that  is,  the  rile  type  abbreviation  suggests  the  type  of  file.  Second, 
they  were  chosen  to  conform  to  the  rile  naming  conventions  in  use  on  the  greatest  percentage 
of  computers;  this  particular  set  can  be  used,  in  conjuncdon  with  the  PC-DOS  and  MS-DOS, 
OS/2,  UNIX,  and  VAX/VMS  operating  systems.  Third,  they  can  be  abbreviated  if  necessary 
for  machines  that  severely  restrict  the  number  of  characters  in  a  file  name;  the  second  and 
third  characters  of  each  designator  form  a  unique  pair.  These  choices,  therefore,  facilitate  the 
portability  and  ease  of  use  of  the  BACKSCAT  code. 


American  National  Standard  Programming  Language  FORTRAN ,  ANSI  X3.9-1978  (New  York, 
NY:  American  National  Standards  Institute,  1978),  Page  12-3. 
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Table  B-l.  BACKSCAT  FUe  Types 


Hie  Structure 

File  Type 

Configuration  Hies 

.CFG 

Lidar  System  Files 

.LDR 

Atmospheric  Structure  Files 

.ATM 

Propagation  Profile  Files 

.PFL 

Resonant  Absorption  Profile  Hie 

.RES 

Internal  Aerosol  Data  Files 

.AER 

Internal  Aerosol  Scaling  Hie 

.SCL 

Simulation  Log  Files 

.LOG 

Simulation  Output  Files 

.DAT 

i 
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APPENDIX  C  Formats  of  Data  Files  Used  by  BACKSCAT 

This  appendix  describes  the  internal  structure  of  files  used  by  the  BACKSCAT  program. 
The  program's  interactive  interface  normally  frees  the  user  from  concerns  about  file  formats; 
thus,  this  information  is  supplied  primarily  to  facilitate  diagnosis  of  corrupted  files. 

C.l  Configuration  Files 

The  configuration  file  contains  all  configuration  information  for  the  simulation.  Bach  con¬ 
figuration  includes  a  default  run  set,  relevant  information  about  the  host  computer  and  operating 
system,  and  the  default  status  of  several  program  options.  Figure  C-l  shows  the  format  of  the 
configuration  file. 


1.  <character> 

1  <character> 

3.  <character> 

4.  <character> 

5.  <ctaracter> 

6.  ccbjracteo 

7.  < logic 'u>  <iog;caI>  <logical> 

8.  < integer >,  <integer>,  <imeger> 

9.  <integer>,  cintegeo  [,<integer>, 

10.  <imeger>.  <integer>  [,<integer>. 


Figure  C-l.  Configuration  File  Format 


name  of  UDAR  system 
name  of  atmospheric  structure 
name  of  propagation  profile 
name  of  molecular  absorption  profile 
name  of  simulation  log 
name  of  simulation  output  data  set 
configuration  flags  for  display 
screen  length,  width,  and  margin 
.)  number  of,  and  character  sequences  for 
ANSI  escape  characters 
.]  number  of,  and  chancier  sequences  for 

ANSI  SET  GRAPHICS  RENDITION  chmeter 


It  is  important  to  note  that  the  data  must  be  input  in  the  specified  older  and  in  the  correct 
format;  however,  the  text  accompanying  the  data  is  optional. 

The  first  six  entries  constitute  the  default  simulation  run  set  these  names  are  convened  to 
the  appropriate  operating  system  file  names  by  the  simulation  by  a  procedure  that  is  dependent 
upon  the  operating  system. 

The  four  LOGICAL  flags  in  the  seventh  line  control  atmosphere  display  (T  for  normal  display 
or  F  for  inverted  display),  propagation  mode  (T  for  atmospheric  model  or  F  for  propagation 
profile),  terminal  type  (T  for  ANSI  terminal  control  or  F  for  a  generic  non-ANSI  terminal),  and 
Rayleigh  scattering  mode  (T  if  Rayleigh  scattering  is  on  or  F  if  Rayleigh  scattering  is  off)- 

The  eighth  line  specifies  the  screen  dimensions  and  the  margins  for  the  status  display  area. 
The  screen  length  must  be  right  justified  in  columns  1-4,  the  screen  width  in  columns  6--9,  and 
the  margins  in  columns  11-14.  Commas  should  be  placed  in  columns  S  and  10. 

The  remaining  two  lines  contain  the  character  sequences  to  be  written  to  the  terminal  for 
the  ANSI  ESCAPE  (line  9)  and  SET  GRAPHICS  RENDITION  (line  10)  characters.  Each  fine 
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mutt  have  the  number  of  characters  in  the  corresponding  string  right  justified  in  columns  1—4; 
the  collating  positions  of  the  characters  are  right  justified  sequentially  beginning  in  columns  6-9 
and  continuing  with  columns  11-14,  16-19,  21-24,  and  26-29  as  required  to  accommodate  a 
maximum  of  five  characters.  Commas  should  be  placed  in  the  intervening  columns,  and  unused 
fields  before  column  30  should  be  blank. 

C2  Udar  System  Files 

Each  lidar  system  file  contains  all  parameters  related  to  a  particular  lidar  system.  The 
parameters  must  appear  at  the  beginning  of  consecutive  records  of  the  file  in  the  specified  order. 
Figure  02  shows  the  format  of  die  lidar  system  file. 

1.  <real>  Wavelength  (job) 

2.  <real>  Transmitted  Pulse  Energy  (J) 

3.  <real>  Pulse  Duration  (joec) 

4.  <real>  Receiver  Aperture  Diameter  (cm) 

3.  <real>  Receiver  Obscuration  Diameter  (cm) 

6.  <real>  Minimum  Range  (Km) 

7.  <real>  Maximum  Range  (km) 

8.  <real>  Range  Resolution  (km) 

9.  <real>  Sensor  Height  (km) 

10.  <real>  Viewing  Azimuth  Angle  (degrees) 

11.  creel  >  Viewing  Elevation  Angle  (degrees) 

Figure  C-2.  Lidar  System  File  Format 

The  comments  describing  each  element,  shown  in  the  example,  are  optional. 

C.3  Atmospheric  Structure  Files 

The  atmospheric  structure  file  contains  the  specifications  required  by  the  AFGL  atmospheric 
model.  Figure  03  shows  the  format  of  atmospheric  structure  files;  data  must  appear  in  the  exact 
sequence  shown,  and  all  text  entries  must  be  entered  un  UPPER  CASE  TYPE. 


1.  <real>  Boundary  Layer  Height  (km) 

2.  <real>  Height  of  Troposphere  (km) 

3.  <real>  Height  ot  Stratosphere  (km) 

4.  creel>  Height  of  Atmosphere  (tan) 

5.  <character>  Boundary  Layer  Aerosol 

6.  <chamcter>  Stratosphere  Extinction  and  Profile 

7.  <chancter>  Upper  Atmosphere  Aerosol  Level 

8.  <chaxacter>  Aerosol  Distribution  Profile 

9.  <real>  Boundary  Layer  Humidity  (%) 

10.  <real>  Troposphere  Humidity  (%) 

11.  <real>  Visibility  (tan) 


Figure  03,  Atmospheric  Structure  File  Format 
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The  first  four  records  in  the  atmospheric  profile  input  file  contain  the  upper  limits  of  each  of 
the  layers-  The  values  in  the  atmospheric  structure  file  supplied  with  the  BACKSCAT  program 
am  correct  for  the  atmospheric  data  supplied  with  the  program;  men  changing  these  values  must 
supply  atmospheric  data  that  is  correct  for  the  modified  values. 

The  boundary  layer  aerosol  on  foe  fifth  record  may  be  URBAN,  RURAL,  MARITIME, 
PCEhNIC,  TRQP0SPHERXC,  ADVECTXVE  FOG,  or  RADIATIVE  FOG. 

The  stratosphere  aerosol  extinction  and  profile  are  defined  on  the  sixth  line.  If  both  foe 
aerosol  and  profile  am  "background",  this  line  should  contain  the  key  word  BACKGROUND.  For 
volcanic  profiles,  columns  1-5  contain  the  extinction  (FRESH,  AGED,  or  BKGD . ,  right  justified, 
columns  6-10  contain  foe  profile  (MOD.,  (moderate)  HIGH,  EXT.,  (extreme)  or  BKGD . ,  alio 
right  justified,  column  11  must  be  blank,  and  columns  12-15  must  contain  the  key  word  vote 
(volcanic). 

The  upper  atmosphere  aerosol  on  line  seven  may  be  either  normal  or  extreme. 

The  eighth  record  of  the  file  contains  the  aerosol  distribution  profile,  which  may  be  either 
SPRING/ SUMMER  or  FALL/WINTER. 

The  boundary  layer  and  troposphere  (relative)  humidity  are  given  in  percent  on  foe  ninth 
and  tenth  records  respectively;  thus,  a  value  of  70.0  corresponds  to  70%  relative  humidity. 

The  last  record  contains  the  visibility  in  kilometers  if  it  is  positive,  or  a  uniform  extinction 
coefficient  if  it  is  negative. 

C.4  Propagation  Profile  Files 

The  propagation  profile  file  may  be  either  read  or  written  by  the  simulation,  depending 
upon  the  simulation  mode.  In  either  case,  it  is  a  table  whose  columns  are:  (1)  altitude,  (2) 
aerosol  extinction  coefficient,  (3)  aerosol  absorption  coefficient,  (4)  aerosol  scatter  coefficient,  (5) 
aerosol  backscatter  coefficient,  (6)  molecular  extinction  coefficient,  and  (7)  molecular  backscatter 
coefficient.  The  table  must  be  arranged  such  that  each  record  after  the  first  ii  at  a  greater  altitude 
than  the  previous  record.  A  maximum  of  100  records  are  permitted. 

C.5  Molecular  Absorption  File 

Hie  user  may  supply  a  two-column  file  to  specify  absorption  by  molecular  resonances  as  a 
function  of  altitude.  The  first  column  of  this  file  must  contain  the  altitude  in  kilometers  and  the 
second  column  must  contain  the  absorption  coefficient  in  inverse  kilometers.  The  data  must  be 
arranged  such  that  the  altitudes  are  in  increasing  order,  that  is,  the  aldtude  on  each  record  must 
be  greater  than  foe  altitude  on  the  previous  record. 

C.6  Simulation  Log  Hies 

The  simulation  generates  a  log  file  when  it  is  executed.  This  tog  file,  written  in  standard  text 
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format,  allows  the  execution  of  the  simulation  to  be  traced.  The  actual  syntax  of  die  simulation 
log  is  dependent  upon  the  simulation  configuration  options;  thus,  the  information  written  to  the 
log  is  documented  in  the  log  itself. 

C.7  Simulation  Output  File 

The  simulation  output  is  a  standard  data  file,  in  text  format,  containing  a  table  of  the  lidar 
propagation  and  backscatter  information  as  a  function  of  range.  The  columns  of  die  table  are:  (1) 
-range,  (2)  altitude,  (3)  cumulative  attenuation,  (4)  lidar  backscatter,  (S)  normalised  backscatter, 
(6)  range-compensated  backscatter,  and  (7)  normalized  range-compensated  backscatter.  This 
data  may  be  plotted  using  most  commercially  available  plotting  programs. 
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